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Knots are closed space curves embedded in three-dimensional (3-d) space. Classification 

Ch ■ 

of knot types depending on their topology is one of central issues in mathematical knot the- 
ory.^ Several invariances, often in the form of polynomial functions, can be used for that 
purpose, calculation of which, however, is generally not an easy task. The computational cost 
h of the knot invariances rapidly increases with the number n of the crossing in the knot dia- 
gram (Fig. 1 and 2), e.g, at least h ~ n for Alexander polynomial, and h increases faster than 

q ■ any polynomial, i.e., NP(non polynomial) -problem for Jones polynomial. 2) 

The p-colorability problem of knots provides a particular type of the invariance (see below 

>• ■ for its definition). 

o : 

\0 , There seems to be no generic method to find its solution, i.e., a colored configuration with 

O ■ 

t-h ■ satisfied manner, called a coloring class, and to count the total number of them. In this note, 

we present a statistical mechanical formulation of the p-colorability problem of knots, which 
provides an algorithm to find the solution. The method also allows one to get some deeper 
j> ! insight into the complexity of the problem from the viewpoint of the constraint satisfaction 

^ ■ problem. 

h : 

Formulation of the problem. The knot topology can be analyzed by using the knot trajec- 
tory projected onto the plane with conserving over- and under- crossing conditions of local 
components. It is called the knot diagram. As symbolically exemplified in Fig. 1(a), it consists 
of n arcs and n crossing points. In p-colorability problem, we attempt to color each arcs using 
p colors under a particular type of the constraint (see Fig. 1(b) ). Let c, e {0 , • • • , p - 1} 
be the color of ?'-th arc. Each crossing point consists of three arcs, at which a constraint is 
defined locally. Suppose that two arcs i and i + 1 are separated by an overcrossing of an arc k 
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Fig. 1. (color online)(a):A schematic example of the knot diagram of figure-8 knot with « = 5. (b):A satisfied 
colored pattern (coloring class) of the knot diagram. (c):A graphical reprsentation of the local satisfiability 
condition (constraint) on each crossing. A crossing and an arc are represented as a clause node (triangle) and a 
variable node (circle), respectively. Each clause is connected with three involved variable nodes, where a node 
corresponding to the overcrossing arc is distinguished from other two. (d): A random graph corresponding to the 
knot (b). The number of edges connected with each variable node is random. 



at z-th cross. If the following equality holds among the colors of these arcs c,-, c i+ i and q, 

mod(c, + a+up) = mod(2c fe , p), (1) 

the z-th cross is satisfied. The total number N(p; %) of coloring classes which satisfy all 
the constraints imposed at each crossings is known as one of the knot invariants, where 9f 
stands for the knot type. The inequality N(p; 7f ) + N(p; %') guarantees the different topology 
7C ± thus, N(p; 7f) is the knot invariant. However, like any other invariants, the reverse 
is not generally valid, i.e., the value of Nip;^) and N{p\'K') are not always different for 
different knot types. 

Mapping onto constraint satisfaction problems. By introducing two types of nodes, we 
obtain a graphical representation of a knot diagram (see Fig. 1(d) ). Variable nodes (circles of 
Fig. 1(c) and (d)) and clause nodes (triangles) correspond to arcs and crossing points, respec- 
tively, which have the following connection properties, 

(i) Associated with each clause node are three edges connected to variable nodes; within z-th 

clause node, z-th , i + 1-th , and k-th variable nodes are involved (Fig. 1(c)). 

(ii) The number m(i) of edges connected to z-th variable node is random under the constraint 
S?=i w(0 = 3n. 

Hence, graphs corresponding to knots have equal number n of clause and variable nodes, 
and are allowed to have random connecting property patterns (by (ii)) with the condition (z) 
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imposed. The satisfiability function /z,(c,, c,-+i, c*)is defined as, 



1 ( mod(c; + c M , p) - mod(2c fc , p)±0) 

hi(c u c M ,c k ) = < 

( mod(c,- + Cj+i , p) - mod(2c*, p) = ) 



(2) 



The degree of the net satisfiability is quantified by the function 



H = } hi(ci,c i+1 ,c k ). 



(3) 



/=i 



Thus a projected knot diagram is mapped onto a random graph, and the problem of com- 
puting N(p; *K) is mapped onto the problem of counting the number of ground states of the 
Hamiltonian (3); a constraint satisfaction problem on random graph. 

Inplement of simulation. We performed replica exchange Monte Carlo simulations^ of the 
Hamiltonian (3) and estimated the ground state entropy, which corresponds to the logarithm 
of N(p; ( K). For computation of the entropy S (fi) at each inverse temperature point (3 = l/T, 
we used following thermodynamic integral equation, 



where Uifi) represents the internal energy. 

First we generated 3-d random sample trajectories with figure-8 knot type by performing 
a standard Langevin-dynamics (LD) simulation using a closed beads-spring model polymers 
(Fig. 2). Random graphs were obtained by projecting the trajectories onto the plane. We then 
carried out MC simulation explained above for each graph and obtained entropy as a func- 
tion of inverse temperature (Fig. 3). To obtain smooth profiles in Fig. 3, simulations were per- 
formed on 40 different temperatures and gaps are interpolated using histogram reweighting 
technique. The number of coloring classes evaluated as an exponentiation of the ground state 
entropy agrees with the known value, i.e., 25 for a figure-8 knot with p = 5 colors. At this 
example, the number of coloring class with satisfied manner is prosperously evaluated. 

Discussion and perspectives The present method is, in principle, applicable to much more 
complex knots than considered here with larger minimal crossing numbers. The results from 
several different coloring number p, and if necessary combined with other kind of the invari- 
ants, would have a good classification ability of knots. 

So far we have only focused on the ground states of the model Hamiltonian eq.(3). Ground 
states are obviously most important in the context of knot theory with a definite meaning 
corresponding to the invariance. It does not, however, necessarily exclude the possible impli- 
cations of the excited states. For instance, we expect that the internal energy of our model at 
finite temperature may contain some useful information on the knot complexity. 
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Fig. 2. (color online) (a):Conformation of figure-8 knot used to perform trial computation. The correspond- 
ing knot diagram has 26 crossings. (b):A realization of the minimal crossing points; 4 of figure-8 knot. The 
diagram(a) with n — 26 and that in Fig.l with n — 5 can be transformed to the diagram (b) with n = «„,„, = 4 

via Reidemeister moves. 
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Fig. 3. (color online)Computed entropy profile for n = 26 system (Fig. 2(a)) as a function of inverse tem- 
perature (circles). Error bars result from the average of 10 independent simulations. For comparison the result 
of minimal crossing number (« = 4) system (Fig. 2(b)) is also shown (crosses), which was obtained by exact 
enumeration of Hamiltonian (3). At low temperature, both converge to the value log 25 within the accuracy of 
the numerical error, see the close up (inset). 



The coloring problem on random graphs often show characteristic phase transition behav- 
ior involving both statics 8 ' and dynamics 9 ' from simple structureless phase to glassy phase 
with many metastable states. 4 ' 7 'lt has been intensively studied from the spin glass perspec- 
tive 5 ' and often discussed in relation to computational complexity. 4 6 ' A unique feature in the 
present system comes from the fact that the ground states of our model are connected to the 
topological invariance. Therefore, we can control the apparent conformational complexity in 
arbitrary ways, i.e., the crossing number n in the knot diagram, while keeping the ground 
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states invariant via Reidemeister moves 1} (See the caption in Fig. 2). The model surely ex- 
hibits extensively separated ground states, but its total number does not grow exponentially 
with the system size n. This fact means that the model does not exhibit clustering and con- 
densation transitions on solution space structure 7 ) in literal terms. Detailed investigations on 
such a model system may be interesting towards better understanding of the glassy properties 
with rugged landscapes 10) and empirical hardness of searching problems. 11_13) 
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